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Spin transport properties at finite electric and magnetic fields are studied by using the generalized 
semiclassical Boltzmann equation. It is found that the spin diffusion equation for non-equilibrium 
spin density and spin currents involves a number of length scales that explicitly depend on the 
electric and magnetic fields. The set of macroscopic equations can be used to address a broad range 
of the spin transport problems in magnetic multilayers as well as in semiconductor heterostructure. 
A specific example of spin injection into semiconductors at arbitrary electric and magnetic fields is 
illustrated. 
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Spin transport phenomena have been identified in various layered systems including magnetic multilayers, magnetic 
tunnel junctions, and semiconductor heterostructure. The partial list of these phenomena includes the spin injection- 
detection in metals jj| , giant magnetoresistance 10 , enhanced Gilbert damping Q , spin angular momentum transfer 
^|, ^| and spin current induced spin excitation jj, ferromagnetic resonant , tunnel magnetoresistance || , and spin 
Hall effect in semiconductors Jk , [n| . Most recently, there are emerging interests in the spin injection and detection 
in semiconductor spintronics [O . Several theoretical models based on a macroscopic drift-diffusion equation for each 
spin channel has been proposed [jl3], |l4|, [ll| . 

The central observation of all these interesting spin-dependent transport phenomena is the key role played by 
the electron spin. Unlike the physics of strongly correlated systems, most of the spin transport phenomena can 
be understood in terms of the semiclassical description of electrons in the diffusive scattering limit as long as the 
spin degree of freedom is properly being taken care of quantum mechanically. In this letter, we construct a set of 
generalized spin dependent macroscopic equations from the spinor form of the Boltzmann equation. We have found 
that these macroscopic equations can self-consistently determine the charge density, spin density, charge current 
and spin currents. The above mentioned spin dependent phenomena can be naturally understood by solving the 
macroscopic equations with various experimental boundary conditions. In particular, one can predict spin transport 
in semiconductor heterostructure at arbitrarily strengths of electric and magnetic fields. 

Let us start with the spinor form of the semiclassical Boltzmann equation that describes the equation of motion for 
spin-half conduction electrons in an external electric field E and an effective magnetic field H e 
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where F is the distribution function that depends on the position r, momentum p (or velocity v = p/m*), and time 
t. H e is the effective magnetic field, H is the spin dependent Hamiltonian which can be written in spin space as 
H = —fiBC ■ H e . The right-hand side of the equation is the collision term. Although the validity of the semiclassical 
Boltzmann equation is limited to the wave packet picture of the electron motion, the quantum mechanical nature of 
the spin has been explicitly taken into account by the spinor form of the distribution function and the spin dynamics 
is included in the commutator [H, F] for the spin variable. 

Since we are interested in spin- dependence of the transport, we will make a few simplifications to Eq. (1). First, 
the effect of the magnetic field on the orbital motion is discarded, i.e., we have set (v x H e ) • V V F = 0; noticeably, 
this orbital term is in fact needed if one is interested in the ferromagnetic resonance under an applied AC and DC 
magnetic fields j8j. Second, we take the collision term in a simplest s-wave form by using momentum-independent 
relaxation time approximations for both non-spin-flip and spin-flip scattering, i.e., 
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where r and r s f are the momentum and spin relaxation times, F = (\/An) j dCl^F^r, v) is the angular average over 
the k-space, I is the 2x2 unit matrix in spin space, and Tr represents the trace over spin space. In ferromagnetic 
materials, the momentum relaxation time r is also spin-dependent; this can be easily generalized by using a 2 x 2 matrix 
for t. For the notation simplicity, we write it as a number. Finally, we should consider only the layered structure so 
that the spatial dependence of the distribution is limited to one-dimension; this is a practical approximation since the 
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three-dimensional inhomogeneity will lead to a complicated spatial-dependence of the current density that makes the 
calculation of the conductance extremely tedious. Let us denote the direction perpendicular to the layered structure 
as ^-direction. 

With above simplifications, Eq. (1) becomes 
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where we have introduced a spin-dependent distribution vector f by separating the distribution function into spin 
dependent and spin independent parts, F = fol + cr ■ f, so that the commutator (l/ih)[H, F] = {^b/K)^ • (H e x f). 

Equation (3) is known as the generalized Bloch equation (ll|. One of major tasks in the spin transport theory is to 
solve this integro-differential equation with desired precision for various physical structures. Here, we do not intend to 
solve exactly the distribution function as a function of the momentum, position and time for a detail structure. Rather 
we construct a set of macroscopic equations from the above equation by integrating out the momentum variable so 
that the spatial dependence of the macroscopic spin dependent transport coefficients can be established for a wide 
range of interesting physical systems. 

The first set of the macroscopic equations are obtained by directly summing over the momentum (velocity). By 
separately writing down the spin dependent and spin independent parts of Eq. (3), we have 
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where we have defined the charge density n n (x,t) = ^v^F, the spin density m(x,t) = ^v^ 1 ^ F), the charge 
current density jo(x, t) = J2 V Tr^F), and the spin current density j m {x,t) = o~v x F). The above equations 

are exact macroscopic equations for charge and spin: Eq. (4) is the law of charge conservation, and Eq.(5) governs 
the motion of spin in the presence of spin relaxation due to spin flip scattering and spin precession due to a magnetic 
field. To determine four macroscopic quantities (no, jo, m and j m ) we need to construct two more equations; this 
can be furnished by timing v x at the both sides of Eq. (3) and summing over v. Again, after separating the spin 
dependent and spin independent parts of Eq. (3), we find 

+ _5npg 1 t) + £^ no(M) = (6) 

djj^ + fagM) + 2- x = 

dt dx m* n t 

In deriving the above closed form, we have used the "mean field" approximation v x ^ « v% ^ v ^ . In another 
words, we only consider the first two polynomials (zeroth and first orders) of the distribution function in the expansion 
respect to the direction of the momentum. This approximation has been used for most of diffusion theories and it 
is valid as long as the system is not too much anisotropic. In principle, however, we should keep higher orders of 
polynomials and one requires to use more macroscopic variables such as Q„ = ^2 v v x Tr(cr F) where n — is for 
the spin (number) density, n = 1 is for the spin current density, n = 2 is for the spin energy density, and so forth. 
Truncating the expansion at n = 2 makes the equations for the macroscopic variables in the closed form, but it 
sacrifices the accuracy in the solution of the Bloch equation for anisotropic systems. However, we argue that inclusion 
of higher orders does not in fact improve much of the accuracy of the calculation since we have already made an 
approximation in deriving the Bloch equation by assuming a k independent relaxation time. 

Equations (4) and (6) describe the non-equilibrium charge density and charge current density. Alone with the 
Poisson equation, Eqs.(4) and (6) completely determine the charge carrier distribution in a layered structure in the 
presence of arbitrary electric fields. A well-known example is the semiconductor p-n junction where the electric field 
in the depletion region is quite large, thus the charge diffusion should be calculated by keeping all orders of the electric 
field in solving Eqs. (4) and (6). Here our focus is on the spin part of the macroscopic equations, Eqs. (5) and (7). 

The spin diffusion equation can be readily derived by expressing the spin current density in terms of the spin density 
from Eq. (7) and by inserting the resulting expression into Eq. (5). Let us first consider the time independent case. As 
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Eq. (7) contains a cross product of the spin current density vector and the magnetic field, it is convenient to separate 
the components parallel and perpendicular to the magnetic field, i.e., 

-xdm z eE x r 



and 
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where we have assumed a uniform magnetic field in the z-direction, and defined the perpendicular components of the 
spin current density — jf n ± ij^ and the spin density m± = m x ± im y . By placing them into Eq. (5), we arrive at 
the spin diffusion equations for the longitudinal and transverse spin density, 
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where S\gn(E x ) = 1(— 1) for E x > (<)0. We have introduced a number of relevant lengths: the usual spin-flip 



scattering length X sc if = \Jv 2 TT s f, the length associated with acceleration of electrons by the electric field Xe = 
m*v 2 /e\E x \, the length due to non-uniform electric field (or due to charge accumulation in the depletion layers) 
A^ 1 = \{\/E x )dE x /dx\, the spin precession length induced by the magnetic field A/j = h^Jv 2 /2fiBH e , and the 



precessional spin density decay length A s = y hTV 2 /2[iBH e . 

The spin diffusive equations, Eq. (10) for the longitudinal spin density and Eq. (11) for the transverse spin density, 
are our central results. For the longitudinal spin density, the spin diffusion is not affected by the magnetic field since 
we have disregarded the electron orbit motion. The electric field as well as its spatial derivative have significant effects 
on the spin diffusion for semiconductor heterostructure where the interface regions are expected to possess large space 
charges and thus large electric fields. For example, for a non-degenerate semiconductor, m*v 2 s» k E T, and thus A^ 
is only about 100 A for a typical electric field of 2.5 x 10 4 (V/cm) at T=300K. In a moderated doped interface, the 
gradient of the electric field in the depletion layer is large so that A c is also of the order of 100 A. Such small Xe and 
A c make the diffusion length highly anisotropic: the spin diffusion is much easier when the direction of the diffusion 
is same as the direction of the electric field; otherwise, the spin diffusion is prohibited. For a metallic layer, however, 
the electric field effect is small since A_e is typically of the order of 10/im if we take m*v 2 « 1 (eV) (a fraction of the 
Fermi energy) and E x is of the order of 10 3 (V/cm) for the extremely high current density. Thus, the commonly used 
spin diffusion equation which neglects the electric field is justified for metallic multilayers. 

We point out that a similar spin diffusion equation for the longitudinal spin density has recently derived from the 
spin-dependent drift-diffusion equations by Yu and Flatte [jl3|, and Albrecht and Smith jl4| with no magnetic fields. 
Our derivation from the Boltzmann equation generalizes their results significantly. First, we are able to address the 
longitudinal and transverse spin diffusion at any magnetic field; the spintronics application is inevitably involving 
the magnetic field. Second, we include the additional length scale A c due to non-uniformity of electric fields at the 
depletion layers; this effect has been neglected all together in the previous treatment. Third, we have started from 
the semiclassical Boltzmann equation so that we have a controlled access on the approximations made in deriving 
the spin diffusion equations; the deviation of the macroscopic equations from the exact solutions of the Boltzmann 
equation can be compared in various cases Jl^j . 

We now turn to a discussion on the transverse spin density and spin current density. The spin diffusion equations 
involve two more length scales associated with the magnetic field in addition to the length scales in the longitudinal 
diffusion equation. The total spin current density decaying length would depend on the detail magnitudes of these 
lengths scales. The magnetic field enters in two places: one is to compete with the spin flip scattering (notice the 
positive sign of 1/A 2 ,) and the other is a precessional decay of spin density due to scattering and precession (notice 
the imaginary i). To give out a quick estimation of the magnetic field on the transverse spin diffusion, we take an 
example of the non-degenerate semiconductor again. We find that the magnetic field of the order of 1000 (Oe) can 
make and X s comparable to X sc if- Therefore, one must consider the spin diffusion by including both finite electric 
and magnetic fields in a realistic spintronics device. 
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Our generalized spin diffusion equations in the presence of finite electric and magnetic fields can be used to a broad 
range of spin transport phenomena. By taking the limiting cases, many previously known spin diffusion phenomena 
are recovered. The giant magnetoresistance for currents perpendicular to the plane of the layers is a simplest example 
where both electric and magnetic fields are assumed small [0. The phenomenon of the spin angular momentum 
transfer with non-collinear magnetization vectors of two ferromagnetic layers can be studied by taking the limit of a 
large internal magnetic field and a vanishingly small electric field fisf| . The ferromagnetic resonant experiments can 
be analyzed via the reduction of the diffusion constant due to a large magnetic field [pi . One can also derive a solution 
of the time-dependence of the spin accumulation in magnetic multilayers by keeping the time dependence diffusion 
term in Eq. (11). 

Finally, we illustrate an important application of our spin diffusion equation in a ferromagnetic/tunnel bar- 
rier/semiconductor heterostructure. To make our calculation much simpler, we assume a uniform electric field near 
the interface of the insulator barrier and the semiconductor; a detailed calculation on the self-consistently determined 
non-uniform electric fields due to depletion potentials and the external bias will be published elsewhere. In the present 
case, the solutions of the longitudinal and the transverse spin density, Eqs. (10) and (11), are simple exponential func- 
tions with the exponential decaying factors given by the above mentioned length scales. To determine absolute values 
of the spin density and spin polarization of the current, we need to use boundary conditions across the insulator 
barrier. Let us introduce a spin-dependent tunnel resistance across the barrier as R = Rq(I + Per ■ Mp) where Mj? is 
the unit vector in the direction of the ferromagnetic magnetization and P is the tunneling polarization. The conditions 
for continuous of spin density and spin current density at the interface between the semiconductor and the insulator 
barrier completely determine the prefactors of the exponential solutions. 
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Fig. 1: Spin-polarization of the current as a function of position in the semiconductor for different electric and magnetic fields. 
The magnetic field is applied at z-axis and the ferromagnet is polarized at y direction for all panels in the figure. P x = jm/jo and 
P y = jm/jo are x and y components of the spin polarization of the current, respectively. The parameters are: the spin-diffusion 
length of the semiconductor is \ s df = 1.0/im, the critical electric field and magnetic field are defined as E c = ksT /e\ a df 
and H c — h/2fiBT s d, the resistance ratio psemi^sdf /Ro = 0.1 (where p S emi is the resistivity of the semiconductor), and the 
polarization of the barrier resistance is P = 0.6. 



In Fig. 1, we show the non-equilibrium spin polarization of the current density as a function of the position in 
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the semiconductor in different electric field and magnetic fields. Here we assume Mj? is at the y direction and the 
magnetic field is applied at the z-direction. At zero magnetic field, the electric field can enhance or reduce the spin 
density, depending on the direction of the field, as pointed out in |H| |l4| . The magnetic field has dramatic effects 
on the spin density and the spin polarization; the spin diffusion length is shorten significantly even for a moderate 
magnetic field, say, 1000 (Oe). On the other hand, one also develops a spin current in z-direction which was absent 
without the magnetic field. This x-component of the spin density and spin current comes from the rotation of the 
injected spin (originally at y direction) by the magnetic field. 

In summary, we have presented a set of macroscopic equations for charge density, spin density, charge current 
and spin current at finite electric and magnetic fields. These equations are the generalization of previous useful spin 
diffusions at various limiting cases. These equations can be applied to a wide range of the diffusive spin transport 
phenomena. This work is supported by National Science Foundation (DMR-0076171). 
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